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Abstract 

With a minimally modified dispersion relation for neutrinos, we reconsider the constraints on 
superluminal neutrino velocities from bremsstrahlung effects in the laboratory frame. Employing 



> 

\^ \ both the direct calculation approach and the virtual Z-boson approach, we obtain the generic decay 
(N 

■ width and energy loss rate of a superluminal neutrino with general energy. The Cohen-Glashow's 



analytical results for neutrinos with a relatively low energy are confirmed in both approaches. We 
employ the survival probability instead of the terminal energy to assess whether a neutrino with a 
given energy is observable or not in the OPERA experiment. Moreover, using our general results we 
perform systematical analyses on the constraints arising from the Super-Kamiokande and IceCube 
experiments. 
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I. INTRODUCTION 



The OPERA collaboration has measured the velocities of muon neutrinos {1/^) with mean 
energy about 17.5 GeV, which travel from CERN to the Gran Sasso. It was found that the 
muon neutrinos travel at a speed faster than light (so called propagating superluminally) , 
and the relative difference of the neutrino speed Vi, with respective to that of light c in the 
vacuum was measured to be 

S'v, = "^^^ = (2.37 ±0.32 (stat.)tSS (sys-)) x 10"' ■ (1) 

Also, the OPERA data indicate no energy dependence [1]. The OPERA result has been 
confirmed by a test performed using a beam with a short-bunch time-structure allowing 
to measure the neutrino flight time at the single interaction level [1]. Interestingly, this is 
compatible with the MINOS results [2] and the earlier short-baseline experiments [3] . From 

the theoretical point of view, many groups have already studied the possible solutions or 
pointed out the challenges to the OPERA anomaly [4-43]. For an early similar study, see 
Ref. [44]. 

The OPERA experiment would strongly imply new physics beyond the traditional special 
relativity if it could be confirmed by the future experiments. If Lorentz symmetry is broken 
hardly [45-47], i.e., there is a preferred frame of reference, we do have two strong constraints. 
The first one comes from bremsstrahlung effects [16]. A superluminal muon neutrino with 
S'v,y given in Eq. (1) would lose energy rapidly via Cherenkov-like processes on their way 
from CERN to the Gran Sasso, and the most important process is — > f^e+e". Thus, the 
OPERA experiment would not be able to observe muon neutrinos with energy in excess of 
about 12.5 GeV [16]. The second one arises from pion decays [17, 20, 22]. The superluminal 
muon neutrinos could not gain energy larger than about 5 GeV from the processes, tt'^ — > 
/i"*"//^ and /i — )■ u^eue [20]. However, these constraints do not apply to other proposals which 
can explain the OPERA anomaly as well [24, 29, 34, 36, 39-42]. 

In Refs. [16, 17, 20, 22], the authors considered the following simple dispersion relation 
for neutrinos as a result of hard Lorentz violation 

El=pJ + ml + 6pJ , (2) 

where E,^ and Pi, are respectively the neutrino energy and momentum, and S ~ 25' v^. For 
simphcity, we shall use S — 5 x 10~^ in the following. The discussions on pion decays are 
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simple, and one can easily confirm the previous results [17, 20, 22] . However, the results for 
bremsstrahlung effects are a bit subtle. Prom private communications with our colleagues 
in October 2011, nobody seems to have succeeded in reproducing the results of Cohen and 
Glashow in the approximation of four- Fermi interactions (see also attempts in Refs. [37, 38]). 
In addition, for neutrinos with much larger energy E > Mz/\/5, their decay width and 
energy loss rate have not been studied as well. 

In this paper, assuming the simple dispersion relation given in Eq. (2), we study in 
detail the electron-position emission process of super luminal neutrinos ly^ — > i/^e"'"e" for 
bremsstrahlung effects in the laboratory frame, and our results can be applied to the other 
processes as well. By a direct calculation, we confirm the Cohen-Glashow's results for the 
low energy neutrinos. We present a general formula for the emission process of superlumi- 
nal neutrinos via a virtual Z-boson exchange, which applies at both high and low neutrino 
energies. Interestingly, we confirm the Cohen-Glashow's low-energy results as well. More- 
over, we systematically analyze the constraints on the neutrino speed from the OPERA, 
Super-Kamiokande and IceCube experiments. 

II. PAIR EMISSION DECAY WIDTH AND ENERGY LOSS RATE 

We study the pair emission process for a superluminal neutrino in two approaches. In the 
first one, we work directly with the low-energy four-Permi interaction that is appropriate 
for momentum transfer much smaller than the Z-boson mass Mz- We recover the results 
reported in Ref. [16]. We also present a new result for the process of on-shell Z-boson 
emission. The latter will be employed in the second approach to work out a formula that 
applies also to neutrinos with general energy. 

A. Direct Calculations 

We start with some kinematical considerations of the pair emission process by a superlu- 
minal neutrino, — >■ i'ij,{£)e~ {pi)e'^ {p2) , where the quantities in the parentheses denote 
the four-momenta 

k^(E,k), e^(E,,e), p,^(Ee-,pi), p2^(E,+ ,p2). (3) 
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The neutrinos are assumed to satisfy the dispersion relation shown in Eq. (2), while the 
electrons fulfil the usual one 



^2 



El = cl\i\' + ctm: 



\pi\^ + ml, 



El 



\p2\'^ + ml, 



(4) 



where nie^u are respectively the electron and neutrino masses, = + 5 is the limiting 
speed of neutrinos in units of the speed of light which we have set to be unity. 

For the energy region of interest here, the lepton masses can be safely neglected, so 
that the dispersion relations are simplified to E = c^\k\, and E^- = \pi\, etc. Noting 
\Pi +P2I < |Pi| + 1^2!; one gets the fraction of the initial neutrino energy carried away by 
the final one as follows 



\k\ 



E^ 
E 



(5) 



where for a given angle 6^, between the moving directions of the initial and final neutrinos 
we have 



/jmax 



COS 9,, 



(6) 



Fig. 1 shows in polar coordinates the relation between i^^^ and 6^, where for illustrative 
purpose we used an unrealistic value 6 = 0.05. The figure implies a factor of 6"^ suppression 
from the neutrino part of the phase space, which will contribute to the decay width. 




FIG. 1: i^^^ is shown in polar coordinates as a function of 9i, for 5 = 0.05. The points (1,0), 
(0,5/2), and (— (^/4, 0) (in rectangular coordinates) correspond to the three kinematical configura- 
tions: k \\ i, k J- i, and —k \\ i, respectively. 

We have done three independent calculations for the pair emission and Z-boson emission 
processes and reached consistent results. In what follows, we present some details about 
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the calculations. The amplitude for the pair emission process at momentum transfer much 
smaller than Mz is 

M = V2GFu(£)^^PLu{k)u{pi)r{{l - 2sI,)Pl - 2s^^Pr)v{p2) , (7) 

where Gp is the Fermi constant, Pl,r = (1 T 7^)/2, spF = sm9w, and u, v are Dirac spinor 
wavefunctions for the particles. By the way, a modified dispersion relation for neutrinos 
may be formulated in Lagrangian field theory in terms of a modified metric tensor, which 
in turn can affect the spin sum of neutrinos and their effective interactions. For a careful 
analysis on such subtleties or model dependencies, see Ref . [43] . 

In this paper, we shall average over the initial neutrino spin states, similar to the Ref. [16]. 
We are aware that the OPERA neutrinos produced from the positively charged pion decay 
have a negative helicity, and the results are roughly the same if one does not average over 
the neutrino spin states [43]. Moreover, our results can be applied directly to the randomly 
polarized neutrinos whose astrophysical sources, for instance, are unknown. In short, the 
spin-summed and -averaged decay width and corresponding energy loss rate are 

r = ^ /^|A4|^rfPS3, (8) 

spins 

<^^lJiK-E),r, (9) 

where dPS^ is the usual three-body phase space measure. Neglecting the neutrino and 
electron masses, which is appropriate for the OPERA experiment, the phase space of the 
electrons can be easily done. We obtain 

r = ^[(l-2.^f + (2.^^f]2^J, (10) 

where the phase space integral J reduces to {q = k — i) 

—^2(k.ie^2k.,i;). (11) 

Using the kinematics defined earlier, we obtain, e.g., k"^ = (5|A;p, and k- i = A\k\\i\. We find 
that it is easier to use the variables £k and y instead of |^| and cosO^, where £k was defined 
in Eq. (5) and y is defined by 2A — {y + y~^)S. The integral becomes elementary 

\k\^ /-i /-^r" 1 

J = / dy / d£, -S{y-' - 1)4 [35A4(1 + £l) - - 2SHl] , (12) 
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where yo — {ci, — 1) / {c,^ + 1) . Working out the integral we get 

mn^cl^^ ^^H^J + l^^^n 1 7 + 20 +20 ij' 

This result is exact with the only approximation being rrie^u = 0. For the energy loss rate, 
one multiplies the integrand of J by —c~^E{l — If.) and obtains 



dx 1927r3c8 
25 

X 



Ig X'j7 3g 1^ '^4 ^3 '^2"^ 



(14) 



^224 

Putting — 1/4 and keeping only the leading term in 5 <^ 1 we recover the Cohen-Glashow 
results [16] 

^ ~ 14^92^ ' ^^^^ 
25 GlEH^ 

— ~ — 5- ■ (16) 

dx 448 1927r3 ^ ^ 

For a supcrluminal neutrino with high enough energy, the process v{k) — )> v{l)Z{p) can 
also take place. Ignoring again the neutrino mass, the spin-summed and -averaged amplitude 
squared is 

^|A^|2 = ./2GF[k-ml + 2k-pl-p\ . (17) 

spins 

The energy-momentum conservation gives 

z'^ - 2zA5-'^ + 1-4 = 0, (18) 

where z — \i\/\k\ is similar to £k in the pair emission process, A was defined in Eq. (6), and 
= {c'uMz) I {E\fb). The above as an equation of z for given cos Qy G [—1, 1] has always two 
real solutions for c^^ > 1, but one of them being larger than unity is non-physical. Requiring 
the other solution to be in the physical region yields the threshold condition, iz < 1, i.e., 
E > c,yMz/VS. Similarly to the case of pair emission, it is simpler to work in phase space 
with z instead of A or cosOi,, 2A = [2; -|- (1 — l\)z~^\5, and the interval for cos^,^ translates 
into 2; e [z- , 2;+] with 

z^ = (2 + 5)5-' [1 - ^1 - (2 + - 4)] , 

z^^l-lz. (19) 



Using the new variables, we have, e.g., 2k ■ £ — c^^£'^(l + z'^ — and 2p ■ k — c^^£^^(l 
z'^ + 1%)5. The decay width and the neutrino energy loss rate are found to be 



167rV2c7 
dE 1 GpE"^^ 



(20) 
(21) 



where the residual phase space integrals are elementary 

h = cl5' r dz [4(1 + z') + (1 - z'f - 24] , (22) 
Ir = clS' r dz [4(1 + z^) + (1 - z^f - 24] {l-z). (23) 



Again the above results only assume me,v — and are exact in parameters £z and 5. Since 
5 <^ 1, we can expand in 5 while holding < 1 fixed and obtain the results in the leading 
order of 5 

r « (4 + 104 - 254 + 114) , (24) 

dE GfE^5^ 



dx 960v^7r 

where now Iz ~ Mz/{E\f6). 



(22 + 354 - I8O4 + 884 + 354) , (25) 



B. Virtual Z Approach 

In the previous subsection we computed the electron-positron emission in the low energy 
limit and the emission of a physical Z boson at sufficiently high energies. In the following 
we present our general results for the pair emission via a virtual Z exchange, ^{k) — > 
p{i)Z*{q) — )> i^{^)f{pi)f{p2), which can be applied at any initial neutrino energy E. One 
can do so purely numerically of course, but we would like to accomplish this in such a way 
that both low and high energy limits can be readily identified. After some manipulations of 
phase space and appropriate reorganization of amplitudes for the subprocesses, we obtain 
upon neglecting the neutrino masses 

dE^ _ 1^ /-SB ^ mSfjrn^) dEi{m^) 
dx~7ri """^^ml-Miy + TlMl dx ' 
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where m* = \f^-, is the usual total decay width of the Z boson, rj(m*) is the decay 
width of the initial superluminal neutrino into a virtual Z boson with an effective mass m*, 
and dEi{m^)/dx is the corresponding energy loss rate. These functions are obtained from 
Eqs. (24) and (25) by replacing Mz with everywhere, i.e., £z ^z* = and 
Gp — )■ G*p = GfM^/ml. Similarly, Tf{m^,) is the decay width of a virtual Z into a pair 
of fermions. The above results can be used to individual channels that are kinematically 
allowed, although in numerical analysis we will show the results summing over all possible 
channels. 

The above Eqs. (26) and (27) serve as a useful function interpolating between the low 
and high energy limits. At energies much lower than Mz but still above the threshold 
for the electron-positron pair, the denominator is approximately equal to M|, and the 
Cohen-Glashow's results in Eqs. (15) and (16) are recovered. In the high energy limit 
E » Mz I V^, the integrals can be worked out in the narrow width approximation (NWA) 
(see, e.g.. Appendix B in Ref. [48]). When summing over all decay channels of the virtual Z 
boson at sufficiently high energy, the results in Eqs. (24) and (25) for the physical Z-boson 
emission are recovered as well. 

III. NUMERICAL ANALYSES FOR NEUTRINO VELOCITIES IN VARIOUS EX- 
PERIMENTS 

A. General Discussion 

We first compare in Fig. 2 numerical results based on Eqs. (26) and (27) with those in 
the low-energy hmit in Ref. [16] and in the high energy hmit using NWA, Eqs. (24) and 
(25). To avoid misunderstanding, we note that while the low-energy result involves only the 
electron-positron emission, the other two include all possible channels that are kinematically 
allowed. Due to the opening of more and more channels one sees a continuous enhancement 
of the decay width (in the left panel) and energy loss rate (in the right panel) . At very high 
energies the NWA result essentially coincides with the full calculation. 

The difference to the low-energy limit can be more clearly seen in the right panel of Fig. 2 
for the energy loss fraction. The mean fractional energy loss is about 0.79 for most energy 
scales until the Z threshold (main peak). At the threshold, the process is dominated by the 
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FIG. 2: Decay width T (in GeV, left panel) and energy loss fraction —{dE/dx)/{ET) (right panel) 
of a superluminal neutrino are shown as a function of the initial energy E (in GeV), using the three 
results: full result (solid-black curve), Cohen-Glashow's result (dotted-blue) and NWA (dashed- 
red). Here (5 = 5 x 10^^, and Eq = 'lm^j\fZ is the threshold for the electron-positron emission. 

production of an on-shell Z which takes away most of the initial neutrino energy. Beyond 
that point, there will be more energy left for the final neutrino. The secondary peak at 
about 30 GeV is due to the setting-in of the mm, rfrf, and ss final states. It is worthy 

noticing that the energy loss rate decreases down to approximately 0.69 in the high energy 
limit where the low-energy approximation is not applicable. 

B. The OPERA Experiment 

Based on the energy loss rate at low energies in Eq. (16) Cohen and Glashow defined a 
terminal energy Et for a superluminal neutrino after travelling a distance of L [16] 

_ 125 G\U^ 

- 448192^^' ^^^^ 
meaning that a neutrino with energy less than Et would not have travelled such a long 
distance. For the OPERA neutrinos with L = 730 km and 6 = 5 x 10"^, they got Et = 
12.5 GeV. We observe from Eq. (15) that a superluminal neutrino with energy Et has a 
decay length, c/F, that is of the same order as the travel distance L. We thus propose to 
use a different statistical measure, the survival probability 

P = e~^^ , (29) 

to decide whether a neutrino is observable or not at the OPERA detector. 
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To get some feeling of numbers, we notice that a superluminal neutrino of mean energy 
~ 17.5 GeV has a hfetime of 1.9 x 10^^ s, which should be compared with its travel time of 
2.4 X 10~^ s if it happens to have not yet decayed during the trip. The survival probability 
can tell in a statistically better way how likely such a neutrino can indeed be observable. 
From the left panel in Fig. 3 one can see that a neutrino has a survival probability of 
0.79, 0.29, 0.09 respectively when it carries a terminal energy 12.5 GeV, mean energy 
17.5 GeV, and an energy of 20 GeV. In particular, a neutrino with an initial energy of order 
Et has a good chance to arrive at the OPERA detector. The right panel of Fig. 3 shows 
that the neutrino lifetime drops drastically as its energy increases. In that case the terminal 
energy serves as an appropriate measure on its observability. 




FIG. 3: Survival probability P (left panel) and lifetime (in seconds, right panel) as a function of 
energy E (in GeV) in the OPERA experiment with 6 = 5 x 10~^ and L = 730 km. 

C. The Super-Kamiokande and IceCube Experiments 

Besides OPERA, the Super-Kamiokande and IceCube experiments also reported obser- 
vation of high energy neutrinos that travel a long distance. The Super-Kamiokande detected 
atmospheric neutrinos that traverse a distance of lO'* km through the Earth (upward-going 
in the detector) over an energy range extending from 1 GeV to 1 TeV [49-51]. The IceCube 
collaboration observed upward-going showers with reconstructed shower energies above 16 
TeV [52], with a baseline length estimated to be at least 500 km. Moreover, it reported the 
upward-going neutrinos from 100 GeV to 400 TeV [53]. For the latter we choose a rough 
travel distance about 10^ km for the neutrinos within a zenith angle of 124° — 180°. 
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We perform an analysis of constraints on the speed deviation 5 of neutrinos from that 
of hght in those two experiments using our general result in Eq. (26). We require that the 
neutrinos detected in Super-Kamiokande and IceCube experiments have travelled a distance 
L that is smaller than n times their decay length Ci^/T at the given energy E, i.e., L < nCj^/T. 
Since T depends roughly on positive powers of E and 5 (e.g., E^5^ at low energy), this gives 
an upper bound on 5 as a function of E. The bound is not sensitive to n as long as n is 
not very different from unity, and a larger n yields a more conservative bound on 5 using 
the same data. Our results assuming n = 10 are shown in Fig. 4 for three baseline lengths 
L — 500, 730, and 10^ km. For instance, for L — 500 km the observations of neutrinos with 
energy in excess of 16 TeV imply 5 < 6.6 x 10^^°, while for L = 10"^ km the observations 
of neutrinos with E > 400 TeV require 6 < 7.8 x 10~^^. These bounds are indeed very 
stringent. 

As both analytical results and Fig. 4 indicate, the bounds on S are not sensitive to the 
baseline length L at a given neutrino energy E but are rather sensitive to E. The simple 
straight-line behavior implies that the extreme high energy regime has not yet been touched 
in those observations. To get some idea of the regime, we plot in Fig. 5 the decay width F as 
a function of E for four values of 6 that correspond respectively to the bounds obtained from 
the experiments OPERA, Super-Kamiokande with {L,E) — (10^ km, 1 TeV), and IceCube 
with {L,E) = (500 km, 16 TeV) and (500 km, 100 TeV). These experiments available on 
the Earth constrain the 5 parameter to such tiny values that it is hard to reach the high 
energy regime with E > Mz / VS. 

IV. CONCLUSIONS 

We revisited the constraints on superhiminal neutrino velocities by assuming a minimally 
modified dispersion relation. We obtained the general decay width and energy loss rate of a 
superluminal neutrino with high or low energy from the bremsstrahlung process via both the 
direct calculation and virtual Z approaches in the laboratory frame. The analytical results 
by Cohen and Glashow were confirmed in the low energy limit. We used a different measure 
to assess whether a neutrino is observable or not in the OPERA experiment. We presented 
new results on the power law for the bremsstrahlung process in the high energy limit. Using 
our general results, we performed systematical analyses on the constraints arising from the 
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FIG. 4: Constraints on (5 as a function of energy E (in GeV) for three baseline lengths L = 
500, 730, and 10^ km from the top to bottom lines. 



10-' 



10-" 



10-= 



10-== 




0.1 10 1000 10' 10' 10» 

E 



FIG. 5: Decay width T (in GeV) as a function of energy E (in GeV) for various 5 = 5 x 10 ^, 
3.6 X 10"^, 6.6 X 10^1°, and 3.0 x lO^^^ from the left to right curves. The horizontal lines indicate 
the decay length equal to the baseline lengths L = 500, 730, and 10^ km from top to bottom. 

Super-Kamiokande and IceCube experiments. 
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